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Abstract. Wc consider the one-dimensional asymmetric simple exclusion process (asep) in 
which particles jump to the right at rate p £ (1/2, 1] and to the left at rate 1 — p, interacting 
by exclusion. In the initial state there is a finite region such that to the left of this region 
all sites are occupied and to the right of it all sites are empty. Under this initial state, the 
hydrodynamical limit of the process converges to the rarefaction fan of the associated Burgers 
equation. In particular suppose that the initial state has first-class particles to the left of the 
origin, second-class particles at sites and 1, and holes to the right of site 1. We show that the 
probability that the two second-class particles eventually collide is (1 -fp)/3p, where a collision 
occurs when one of the particles attempts to jump over the other. This also corresponds to the 
probability that two ASEP processes, started from appropriate initial states and coupled using 
the so-called "basic coupling" , eventually reach the same state. We give various other results 
about the behaviour of second-class particles in the ASEP. In the totally asymmetric case (p = 1) 
we explain a further representation in terms of a multi-type particle system, and also use the 
collision result to derive the probability of coexistence of both clusters in a two-type version of 
the corner growth model. 



1. Introduction 

In the one-dimensional asymmetric simple exclusion process (asep), particles perform continuous- 
time random walks on Z with rate p € (1/2, 1] of jumping to the right and q = 1 — p oi jumping 
to the left. Particles interact by exclusion; attempted jumps to occupied sites are suppressed. 
The resulting process r]t is Markov in the state space {O, •}^, where for a site a; € Z, r]t{x) = • 
indicates that the site x is occupied at time t, and %(a;) = O indicates that it is empty. If p = 1 
then all jumps are to the right and the process is known as the totally asymmetric simple exclusion 
process (tasep). 

A convenient way to construct the process is using a family of independent Poisson processes 
{N{x^ X ± 1), .T e Z}, where for each x, N{x, x + 1) has rate p and N{x, x — 1) has rate q. At the 
time of a point in the process N{x,y), if there is a particle at site x, it attempts to jump to site 
y; the jump succeeds if site y is empty. Two versions of the process started from different initial 
conditions can be jointly constructed under the basic coupling by using the same N{x, y) for both 
versions. 

A second-class particle, denoted ®, is a particle that interacts with particles like a hole, and 
interacts with holes like a particle. That is, a second-class particle jumps right at rate p if there is 
a hole on its right, and jumps left at rate q if there is a hole on its left; if the site on its left contains 
a particle, they exchange positions at rate p, and if the site on its right contains a particle, they 
exchange positions at rate q. 

A second-class particle can be regarded as a discrepancy between two systems evolving together 
under the basic coupling. If the systems start from two initial conditions which coincide except 
at a single site, where one has a particle and the other has a hole, then there will be exactly one 
discrepancy between the two systems for all time, and this discrepancy moves as a second-class 
particle. 
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Consider an initial state in which every negative site is occupied by a particle, every positive site 
has a hole, and there is a single second-class particle at the origin. Let X{t) be the position of the 
second-class particle at time t. For the tasep, Ferrari and Kipnis proved that X(t)/t converges 
in distribution as i ^ oo to a uniform random variable on the interval [—1,1]. The extension of 
the method of [2] to the ASEP is straightforward, but for completeness we give the result and its 
proof (indeed it can be extended much further, to more general asymmetric exclusion processes 
and to a larger class of initial distributions). For the particular case of the tasep, almost sure 
convergence has been proved by Mountford and Guiol [5|, Ferrari and Pimentel [SJ, and Ferrari, 
Martin and Pimentel [S]- 

We then consider systems with two second-class particles. We will be particularly interested in 
the event that one of the second-class particles attempts a jump onto the site of the other one. 
Various conventions are possible for what happens when such a jump is attempted; the particles 
may annihilate, or coalesce, or one may have priority over the other, or no change at all may 
occur. (See also Section [9] for discussion of another relevant possibility). The difference between 
these conventions will not be important for our purposes, since we concentrate on the probability 
that such a jump attempt happens at all. For convenience we will refer to such an event as the 
collision of the two second-class particles. 

We start the process from a deterministic configuration, denoted •®®6, that has all negative 
sites occupied by (first-class) particles, the origin and site 1 occupied by second-class particles, and 
all sites to the right of site 1 empty (in this notation the origin is underlined and • and O denote 
the semi-infinite sequences . . . ••• and OOO ■ ■ • respectively) . We show that the probability 
that the two second-class particles collide is If instead the particles start at distance two 

apart, from the configuration •®0®6, then the corresponding probability is ^^^2 ■ 

These probabilities can also be interpreted as "coupling probabilities" for systems with only one 
type of particle. If two one-type ASEPs are started from the initial conditions ••OO and •0*0 
and evolve together under the basic coupling, then the probability that they eventually reach the 
same state is ^^p'- '^^^ second probability ^^^f applies to the case where the initial conditions 
of the two systems are ••OOO and •OO^O. 

In the particular case of the tasep, where jumps only occur in one direction, various interesting 
alternative representations are possible. Consider a multi-type system in which the initial condition 
has a particle labelled i at site i, for all i S Z. Lower- numbered particles have priority over higher- 
numbered particles; that is, the particles start out in reverse order of priority and every jump 
consists of a particle labelled i exchanging places with a particle labelled j to its right, for some 
i < j. Whenever the particle at site x has priority over the particle at site a; -I- 1, such an exchange 
between sites x and a; -I- 1 occurs at rate 1. From the coupling/collision probabilities described 
above, one can deduce that the probability that at some point particle i overtakes particle i -I- 1 is 
2/3, while the probability that particle i overtakes both particle i -I- 1 and particle i -I- 2 is 1/2. 

The TASEP with two second-class particles can also be mapped onto a two-type version of 
the "corner growth model" associated to directed last-passage percolation. This model has also 
recently been studied by Coupler and Heinrich [T] . The event that the two particles stay separated 
for ever in the TASEP corresponds to the event that both clusters in the growth process grow 
unboundedly, without one surrounding the other. In ^ it is proved that this event has positive 
probability. Using the correspondence with the TASEP, we obtain that the probability is in fact 
1/3; this gives perhaps the first example of a model where such a "coexistence probability" can 
be precisely calculated. 

In Section[2]we define notation and state our main results for the general ASEP, along with some 
discussion. In Section[3]we discuss the particular case of the tasep and explain the representation 
as a multi-type particle system, and in Section [4] we describe the two-type competition growth 
model. Section |5] collects together some known results on hydrodynamics and couplings. Proofs 
of our results are given in Sections |6] [7] and [8| 



COLLISION PROBABILITIES IN THE ASEP 



3 



2. Statement of Results 

Let p e (1/2, 1] and write q = 1—p. We define the two-type ASEP as follows. For i,j e {O, ®, •}, 
we say that i has priority over j and write i ^ j if i — G {O,®} or if i = ®,j = O. Let 
S — {O, ®, 9}^ and consider the generator given on cylinder functions / : S ^ M by: 

(2.1) Lfif^) - J2 {pMni^Uni^+m [/('?"•"+') " fiv)] + [/(?7"^"-') - fiv)]}, 
where 

{r]{z), iiz^x,y 
rj{x), ifz^y 

If the initial configuration rj belongs to {0,#}^, we recover the ASEP with a single type of 
particle. If p = 1 the attempted jumps are only to the right, and we have the tasep. 

The existence of this process can be proven using Liggett 's construction [7] or using the Harris 
graphical construction as in [2111]. See also the comments in Section [s] about the basic coupling. 

Our results concern systems in which all sites to the left of a certain point contain first-class 
particles, and all sites to the right of a certain point contain holes. (This property is preserved 
by the dynamics). As mentioned in the introduction, we use a compact notation for such states: 
• = . . . ••• denotes a semi-infinite string of first-class particles, 6 = OOO . . . denotes a semi- 
infinite string of holes, and an underline indicates the position of the origin. For example, •®®0 
denotes the configuration defined by 

( •, if a; < 0; 
(•®®6)(a;) := < ®, if a; G {0,1}; 

[ O, if a; > 1 . 

(In fact the absolute position of the origin is rarely important for the results, and sometimes we 
may omit the underline, but it is often convenient to fix it). 

The following result extends to the ASEP a result proved for the tasep by Ferrari and Kipnis 

m- 

Theorem 2.1. Consider the ASEP with p e (5, 1], starting from the state •®6. Let X{t) be the 
position of the second-class particle at time t. Then X(t)/t Up in distribution as t ^ go, where 
Up is uniformly distributed on the interval [~(p— q),p~ q]. 

We note that in fact the same method of proof can be used to prove a much more general 
result. Already in [2j, the initial condition can be taken to be the product measure with density 
A on negative sites and density p on positive sites, with X > p. In addition one can extend to 
more general asymmetric exclusion processes, in which the particles (still interacting via exclusion) 
carry out more general continuous-time random walks with some drift 7 > 0. The convergence is 
then to a uniform distribution on the interval [7(1 — 2A),7(1 — 2p)]. 

In the particular case of the tasep, the convergence has been shown to hold almost surely [8], 

We now turn to systems with a pair of second-class particles. The following theorem again 
extends to the ASEP results proved for the tasep in [2\. The proofs are analogous; we include 
them for completeness. 

Theorem 2.2. Consider the ASEP with p g (i, 1], starting from •®®6. Let r be the first time 
that a jump is attempted by one of the second-class particles onto the site of the other (which may 
be infinity). For t < t, let X(t) and Y(t) be the positions of the two second-class particles at time 
t, with X{t) < Y{t). 

(a) For any r e [-{p - q), {p - q)], 

(2.2) lim P(r > t, X{t) < rt < Y{t)) - " " ""^ 



f^oo ^ ' 4p(p — q) 
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(b) 

(2.3) lim -E\(Y{t) - X{t))l{T > t) 



{P - if 
3p 



In display (3.7) of |2], the limit (2.2) was stated for the case p — I, but it has the extra term 
"—ru{r, 1)" which does not belong there. 
The following is our main result: 

Theorem 2.3. Consider an ASEP with p e (^,1]; containing two second-class particles. Let t 
be the first time that a jump is attempted by one of the second-class particles onto the site of the 
other (which may be infinity). 

(a) // the initial condition is •®®0, then 
(2.4) P{t < c^) = 

3p 

(b) // instead the initial condition is •®0®0, then 

l + 2p2 



(2.5) P(t < oo) 



In the TASEP {p = 1) these probabilities are 2/3 in part (a) (an upper bound of 3/4 was obtained 
by |2] in this case) and 1/2 in part (b). The value in part (b) is obtained from that in part (a) using 
a short argument based on conditioning on the first jump and some symmetries of the process 
under this initial configuration. 

These probabilities also represents a coupling probability in the context of the ASEP with only 
one class of particle. Consider two processes started from the states ••o6 and •0#6, and 
coupled using the basic coupling, under which jumps are attempted at the same sites at the same 
times in both processes (see Section [s] for details). At time there are discrepancies between the 
two processes at sites and 1; these discrepancies move around, each behaving as a second-class 
particles, and at time r, the first time that one jumps onto the site of the other, the two processes 
reach the same state and stay coupled thereafter. Hence the probability that they ever reach the 
same state is or 2/3 for the case of the tasep. If instead the initial conditions are ••006 
and •OO*® then the two discrepancies start at sites and 2, and the coupling probability is 
^g^l* , giving 1/2 for the case of the tasep. 

Our proof of Theorem |2 . 3| involves estimating the expected number of particles which are to the 
right of the second-class particle at time t in a process started from the initial condition •®6. As 
a corollary of these methods, we can deduce the following result, which relates to the convergence 
of the environment around a second-class particle: 

Corollary 2.4. Let rit,t>Q be the ASEP with p G (|, 1] started from the state •®6. Let X{t) be 
the position of the second-class particle at time t. Then 

(2.6) Wu, Ei^t{X{t)-l) = \-^~'^ 



and 



1 , P 



(2.7) lim Er^t{X{t) + l, o ■ a ' 

t— +00 ^ ZD 

3. MULTI-CLASS TASEP 

When p — I, the system is "totally asymmetric" and various interesting connections can be 
drawn. 

We consider a system with infinitely many classes of particle. Initially, site i contains a particle 
labelled i, for all z G Z. A particle with label i has priority over a particle with label j ii i < j. If 
the particle at site x has priority over the particle at site a; -I- 1, the two exchange places with rate 
1. Since the process is totally asymmetric, once two particles have exchanged places, they can 



COLLISION PROBABILITIES IN THE ASEP 



5 



never swap back again, so the higher-priority particle stays ahead of the lower-priority particle 
thereafter. 

Let Xi{t) be the position at time t of the particle with label i (so Xi{0) = i for all i). Consider 
i = for example. Particle has lower priority than all the particles starting to its left, and higher 
priority than all the particles starting to its right. Hence Xo{t) behaves just as the position X{t) 
of the second-class particle starting from the configuration #©6. The same is true for any other i 
(shifted by a constant distance i to the right). In particular, from the strong law of large numbers 



mentioned after Theeorem 2.1 we know that the limits 

(3.1) D,^limX,{t)/t 

exist almost surely for all i, and all have the uniform distribution on [—1, 1]. 

Now consider the particles labelled and 1. Again they see stronger particles to their left and 
weaker particles to their right, so their behaviour corresponds to that of the pair of second-class 



particles starting from the configuration •®®6. In particular, we know from Theorem 2.3 a.) 
that the probability of a jump attempt between these two particles is 2/3. Since the process is 
totally asymmetric, the only way for this to happen is that particle attempts to jump right onto 
particle 1 at some point. In this case the two particles exchange places. The same applies for any 
pair of particles i and i -t- 1 . So we obtain the following result giving the probability that a particle 
eventually overtakes its original right-hand neighbour: 

Corollary 3.1. For any i, P{X,{t) > Xi+i{t) for some t) = 2/3. 



From Theorem 2.2 ^a), we can obtain further information about the joint distibution of the 



"asymptotic speeds" Di and -D^+i of neighbouring particles i and i + 1, as defined at (3.11. We 
get that 

P(A < r, A+i >r)= ~ ""^^^ + ''^ = P(A < r)P(A+i > r) 

for all r G [—1, 1]. This expression would be consistent with Di and -Di+i being independent; how- 
ever, since we know that particle i overtakes particle i+l with probability 2/3, such independence 
does not in fact hold. 

From the second part of Theorem |2.3| we can make the following deduction about the probability 
that a particle overtakes both of the two particles starting just to its right: 

Corollary 3.2. For any i, P{X^{t) > ma.x{Xi+i{t),Xi+2{t)} for some t) = 1/2. 

To deduce Corollary |3.2| from Theorem |2.2[b), one can compare the multi-type tasep of this 



section with the two-type tasep of Theorem 2.2 ^b), started from the initial condition •®0®6, 
as shown in Figure [l] 



■2 -10 



Figure 1. Aligned initial conditions for the multi-type tasep (above) and the 
two-type tasep (below) 

Let the two processes evolve together under the basic coupling; potential jumps occur at the same 
site at the same times in the two processes. 

FolUowing the joint evolution of the two processes, one can observe that the site occupied by 
particle (in the multi-type process) always corresponds to the site of a second-class particle (in 
the two- type process). In addition, the other second-class particle is at the site either of particle 1 
or of particle 2, whichever is further to the right at the time. So the first time that one second-class 
particle tries to jump onto the site of the other is precisely the time that particle jumps onto 
the site of whichever is further right of particle 1 and particle 2, i.e. the first moment that particle 
has overtaken both particle 1 and particle 2. 
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A similar argument can be used to show the following; the probability that particle i overtakes 
all of particles i + l,i + 2,...,i + min the multi-type process is the same as the probability that 
the two second-class particles colhde in a two-type process starting from the state #©0 . . . 0®6, 
where there is a string of m — 1 holes between the two second-class particles. We have seen that 
this probability is 2/3 for to = 1 and 1/2 for to = 2. Simulations suggest that the probability is 
in fact 2/(to -I- 2) for all to. 



4. TASEP AND COMPETITION GROWTH 

In this section we describe the two-type growth model which corresponds to the tasep with 
two second-class particles. 

The connection between the tasep and the corner growth model dates back to Rost [TU]. In 
the simplest case, the growth takes place in the positive integer quadrant (where we write 
N = {1,2,...}). At time all the sites of this quadrant are empty. A site {x,y) may become 
occupied if both sites (x — 1,?;) and {x,y — 1) are occupied; once this condition holds, the time 
until site (x, y) becomes occupied is an exponential random variable with rate 1 (independently 
for each site). As a boundary condition one assumes that the sites {Q,x) and (x, 0) are always 
occupied, for each a; e N. 

We will consider an extension of this model, where each occupied site has one of two colours, 
red and blue. A site takes its colour at the moment it is occupied and never changes thereafter. 
In the initial condition, the sites (1, 1), (1, 2) and (2, 1), and all other sites are unoccupied. (The 
colour assigned to the three already-occupied points is not important). 

When a site becomes occupied, its colour is chosen as follows: 

• Any site {l,x) or (x, 1), for a: > 1, becomes blue. 

• The site (2, 2) becomes red. 

• A site {x,y), x > 1, y > 1, {x,y) ^ (2,2), takes the colour of whichever of its parents 
became occupied more recently. 

(Since the delay times for occupation of each site are independent continuous random variables, 
with probability one all times of occupation are different and so there are no ties to break in the 
third case above). 

There are three possibilities for the first event of the growth process: (1,3) becomes occupied 
and blue, or (3, 1) becomes occupied and blue, or (2,2) becomes occupied and red. Each of these 
occurs with probability 1/3. 



1 



12 3 4 



Figure 2 . To the left an example showing a growth process after its first three 
jumps. We can see that G{2,2) < G{3,1), since site (3,2) took its colour from 
(3,1) rather than (2,2). Also clearly G(2,2) < 0(1,3) (since (1,3) is not yet 
occupied). As a result site (2,3) will take its colour from (1,3) and so will also 
be blue; the red cluster will be surrounded and unable to grow further. To the 
right, an example showing a growth process after its first nine jumps. Here one 
can deduce, for example, that the first site to be occupied was (1,3). Given that 
(1,3) is the first to be occupied, the conditional probability that the red cluster 
survives is 1/2 (see the remark at the end of this section). 
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Note that there will certainly be infinitely many points which take the colour blue (for example, 
all the points on the axes (x, 1) and (l,a;)). However, there may or may not be infinitely many 
red points. 

For example, suppose that the first event in the process is that (2, 2) is occupied (red). In this 
case, both (1,3) and (3,1) are occupied later than (2,2). As a result, (2,3) will take its colour 
from (1,3) (its more recently occupied parent) and similarly (3,2) will take its colour from (3, 1). 
So both (3, 1) and (1,3) become blue, and the red cluster has been "surrounded"; there will be no 
further red sites. 

If (2, 2) is not the first site to become occupied, the red cluster may nevertheless become 
surrounded later; alternatively, it may grow unboundedly See Figures [2] and [7] for some examples. 

Theorem 4.1. The probability that the red cluster grows unboundedly is 1/3. 




Figure 3. Two simulations of the growth process up to time 300. The left 
simulation is a rare example where the red cluster grows for a long time (here 
up until time 271.8) but then becomes surrounded. On the right, the red cluster 
looks likely to survive for ever. 



Coupler and Heinrich [IJ recently showed that the probability in Theorem 4.1 was strictly 
positive, without using any tasep representation, by proving that, in the associated directed last- 
passage percolation model, there exists with positive probability an infinite geodesic starting from 
the point (1, 1) and passing through the point (2, 2). 
Theorem 



4.1 



is again deduced from Theore m|2.3[ a). We give the details of the proof in Section 
[9] The main element in the proof of Theorem |4.1| is the representation introduced by Ferrari and 
Pimentel [6J in which a site containing a second-class particle particle ® is replaced by a pair of 
sites in the state [0,#]. Ferrari and Pimentel treated systems with just one second-class particle, 
so that the two-type system with initial condition •®0 was replaced by a one-type system with 
initial condition •0#0. In Rost's correspondence, this equates to a growth process starting from 
the state where site (1, 1) alone is occupied. 



Analogously, the initial condition •®®C) from Theorem 2.3 'a) corresponds to the initial con- 



dition •0#0*C), or to a growth model starting with sites (1, 1), (2, 1) and (1, 2) occupied. 

This hole-particle representation of two second-classs particles survives only until the moment 
when one second-class particle attempts to jump onto another; but this is enough for our purposes, 
since we will be concerned only with the probability of such a collision event occurring. 

In [6], the path of one second-class particle corresponded to the interface between two clusters. 
Here we have two interfaces, between the red cluster in the middle and the blue clusters on either 
side of it. The event that the two second-class particles meet will correspond to the event the two 
interfaces meet. If this happens, the red cluster between the two interfaces is surrounded and is 
unable to grow any further. This will happen with probability 2/3, as in Theorem 2.3 a). See 
Section [9] for full details. 
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Remark. Note that the first site to be occupied in the growth process is equally likely to be (1,3), 
(2, 2) or (3, 1), with probability 1/3 each. As observed above, if site (2, 2) is the first to be occupied, 
the red cluster cannot survive. We know that the overall probability that the red cluster survives 
is 1/3. Hence, given that the first site to be occupied is not (2, 2), the conditional probability that 
the red cluster survives is 1/2. This statement corresponds to part (b) of Theorem 2.3 in the case 
p — 1, (while Theorem 4.1 corresponds to part (a) of Theorem |2.3[) 



By symmetry, it does not 
matter whether (3, 1) or (1, 3) is occupied first. To emphasise the connection, note that if we start 
from the configuration where (1,3) has been occupied in addition to (1,1), (1,2) and (2,1), this 
corresponds to starting the tasep from the configuration •0#00#6, which, under the hole- 
particle representatio n of the second-class particle, corresponds in turn to the initial condition 
9®0®0 of Theorem 2.3 b). Again, see Section [9] for details of the correspondence between states 
of the TASEP and states of the growth model. 



5. Hydrodynamics and coupling 

In this section we state some known results on hydrodynamics and on coupling that will be 
used in the proofs. 

Hydrodynamics. Define the product measures i^A, < A < 1 under which every site is occupied 
by a particle with probability A and a hole with probability 1 — A, independently for different 
sites. Any translation-invariant stationary distribution of the one-type ASEP is a mixture of these 
product measures [7j. 

Suppose the initial distribution puts probability 1 on configurations with asymptotic density 
1 to the left and to the right. (In particular this is the case whenever all sites to the left of 
some point contain particles and all sites to the right of some point contain holes, i.e. the initial 
configuration has the form • a O for some finite string a) . 

Then the process converges to local equilibrium: 

(5.1) lim S{tN)eyrN\f ^ Vu(r,t)f 

for any cylinder function /, where 0^ is the translation by x, S{t) is the semigroup of the ASEP 
with generator given by ( |2.1[ ), [-J is the integer part and for t > 0, r £ M, u{r,t) is the entropy 
solution of the inviscid Burgers equation 

f dtu{r,t) + {p-q)dr(u{r,t){l-u{r,t))) =0 
\ u{r,Q) = I{r < 0) 

which is explicitly given by 

r 1, r<-{p-q)t; 

(5.2) uir,t)^l fc^, -{p-q)t<r<ip-q)t; 

[ 0, r> {p-q)t. 

Basic Coupling. The Harris construction of the process makes use of a family of independent 
Poisson processes ui = {N{x, a; ± 1), x G Z), where N{x, a; -I- 1) has rate p and N{x, a; — 1) has rate 
q . These live in a probability space {^l,T,P). For the process with first-class and second-class 
particles, we write S = {0,®,#}^ as above, and we define a function (f> : il x S x _R+ S 
inductively in time as follows (informally): fix an initial configuration 77 and call rjt = ^{uj,ri,t). 
Then 

(?7t_)^^y, if t £ N{x, y) and i-j{x) < 7j{y) 
rjt- , otherwise 

that is, the particles attempt to jump from a; to ?/ at the epochs of the Poisson process N{x, y) and 
do so if their classes admit. Of course this definition is incomplete, as the epochs of the Poisson 
process are not well ordered, but a simple argument shows that for all t and almost all realizations 
cj, Z is partitioned into finite intervals that do not interact with each other in the interval [0,t\ 



Vt 



(see t3j and [4j for details). The process so defined is Markov and has generator L given in (2.1) 
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The basic coupling between n versions of the process with initial configurations 77 = (ry^, . . . , 77") G 
S" is the process which uses the same oj for all marginals: 



(5.3) 



That is, in the basic coupling the particles attemp to jump from x to y aX the same epochs in each 
marginal, always respecting the classes (at the marginal). 



6. Speed of a second-class particle 



Proof of Theorem 



2.1 



As in the proof of Theorem 1 of |2] , the main ingredient is a coupling 
argument and the convergence to local equilibrium (5.1 1. One wants to show that for re [— (p — 

p — q — r 



(6-1) limP(X(.)>H) 

We will consider ASEP processes rjt and Q started from the initial states rjo = ••(!). and Co = •06. 

Let Jlirjo) be the number of particles of r]t which are to the right of rt, and similarly for J[((^o)- 
We will compute E{J^{i]q)) — E{Jl{(^o)) in two different ways. 

The initial states ryo and Co differ only at the origin, where 770 has an extra particle. If we couple 
the processes using the basic coupling, then there is exactly one discrepancy at all times, and it 
behaves as the single second-class particle of a process started in the state 9®0. Hence 

(6.2) EiJ^im)) - E{J: (Co)) = P{X{t) > rt). 

On the other hand, by shifting the origin by 1, we can see that the number of particles to the 
right of rt in Ct has the same distribution as the number of particles to the right of rt + 1 in rjt. 
Hence 



(6.3) 



E{J^{^o))-E{J^{Co)) = P{vt{\rt^) 



1 



which, by the convergence to local equilibrium in (5.11, is ^ '' 



give (6.1 1 as desired. 



2(p-g) 



SO that (6.2 1 and (6.3 1 together 

□ 



7. Distance between two second-class particles 
In this section we give the proof of Theorem |2.2[ For r G [—1,1] let 

(7.1) Jl-{v) ■■= J2 Vt{x)+Vt{[rt\){[rt\ + 1 - rt) 

x>rt 

be the number of particles at the right of rt at time t; for convenience, we consider that if there 
is a particle at [rt\ a fraction ([rtj + 1 — rt) of it is to the right of rt. 

Proof of Theorem |2.2[ Consider an ASEP rit,t>0 started from the initial condition 7]a = t) = 

••6 = ••06. 

From fi, the first jump of the process is to the state C)0#6 and this jump occurs at rate p. 

Hence using the definition of in (7.1 ) and the Kolmogorov backward equation, we have that 
EJtifl) is continuous and is differentiable everywhere except possibly at points t where rt G Z, 
with derivative 

(7.2) |i?(J[(77) =p[i?(j[ (•0«6) -i?(J[ (••06))] -rE{i^,{lrt\)). 

We couple two versions of the ASEP starting from •0*6 and from ••06, aligned as follows: 

(7.3) •^•^ 
^ ' ••00 

At time zero there are two discrepancies, located at sites and 1. If a jump is attempted from the 
site of one the discrepancies into the other, the discrepancies disappear and the two configurations 
coincide from that instant on. Until that time, denoted by r, the discrepancies behave as the 
positions of two second-class particles, X{t) and Y{t). The difference of the number of particles 
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to the right of rt for the two coupled processes can be expressed in terms of the positions of X{t) 
and Y{t), as follows: 



(7.4) 



J[(«o«6)- J[(««o6) 



1, if T X and X{t) < lrt\ < Y{t) 

rt - [rt\ , if T > t and X{t) = [rt\ < Y{t) 

l-rt+ [rt\ , if T > t and X{t) < [rt\ = Y{t) 

^0, otherwise. 



Inserting (7.4) in (7.2|, 



(7.5) -E{J[m=p[p{X{t) < [rt\ < Y{t)) + {rt - lrt\)P{X{t) = [rt\ < Y{t)) 

+ (l~rt+lrt\)P{X{t) < lrt\ ^Y{t))] - rE{r^tilrt\)) 



Now we compute (7.2) in a different way by associating to J^^tiq) the following martingale: 

Mt = JtiVa) - / {{lrs\ + 1 - rs)wirs\-iiVs) + (rs- lrs\)wirs\iVs) - '"»7s(ksJ)} ds, 
Jo 

where, for a fixed configuration 77, Wxirj) denotes the instantaneous current through the bond 
[x, X + I], namely: 

Wxiv) = P?y(a;)(l - vi^ + 1)) - Qfli^ + 1)(1 - Vi^))- 
This expression together with the convergence to local equilibrium ( |5.1| imply 



(7.6) 



lim -E{Jl{7))) = {p- q)u{r, 1)(1 - u{r, 1)) - ru{r, 1). 
t—i-oo dt 



The same convergence to local equilibrium also gives that 

(7.7) E{f]t{[rt\))^u{r,l) ast^(x. 
Note that 

(7.8) lim P(X{t) = [rt\) = lim P(Y{t) = [rt\) = 0. 

(For example, X(t) represents the discrepancy between two coupled processes starting from the 
states •06 and ••6; the probability of finding this discrepancy precisely at the site [rt\ at time 
t goes to as t ^ 00 since, again by the convergence to local equilibrium in (5.1 ), the probability 
of site [rij being occupied at time t converges to the same limit u(r, 1) for both processes.) 
Putting equations (7.5l-(7.8l together gives 



lim P{X{t) <rt < Y{t)) 



p-q 

P 



u{r,l)[l~{u{r, 1))]. 



The form of the function u from (5.2) then gives (2.2). 



The fact that X(t) < Y{t) whenever t > t, and the relation 



E 



{Y{t) - X{t))^I{T >t)\=Y^ P(x{t) < y < Y{t),r > t) , 



then give 



lim -E 

t—>oo t 



{Y{t) - X{t))l{T > t) 



1 

4p(p - q) 



p-q 



ip-g) 



[ip-q? 



dr. 



giving (2.3) as required. 



□ 



COLLISION PROBABILITIES IN THE ASEP 



11 



Collision probabilities 



Proof of Theorem 2.3 'a). Denote by •®OC) the configuration that has a second-class particle 



at the origin, while the negative sites are occupied by first-class particles and the rest is empty. 
Denote by ••®6 the configuration that has a second-class particle at site 1, while the non- 
positive sites are occupied by first-class particles and the rest is empty. We couple both processes 
and prove later that the positions of the discrepancies initially at sites and 1 behave as X{t) and 
Y{t), the two second-class particles of the theorem, with X{0) = and Y{0) = 1. 

For a configuration 77 with a unique second-class particle and a finite number of first-class 

particles to the right of it, denote by X{t, rf) the position of the second-class particle at time t and 
(1') 

by Jl ijf) the number of first-class particles to the right of X(i,7]) at time t: 

t(2) 



where rjt is the configuration at time t for the process with initial configuration r;o = ?/. 

Define the configuration fj = •®0. For 770 = fj, the Kolmogorov backwards equation gives 



(8.1) 



dt 



In (8.1 1 %®%C) is attained when the rightmost first-class particle at site —1 jumps to 0, the site 



occupied by the second-class particle, and •0®0 arises when the second-class particle jumps 
from to 1, site occupied by the leftmost hole. The law of (jf) does not depend on the actual 
location of X{Q,rj) but on the relative positions of the other particles with respect to ^(0,77). 

So, we are free to change the origin. In particular E (^j['^\%®%0)^ = e(^J^'^\%®%0)^ and 
2i;(jf ^(•®6)) = £^(jf' (•®06)) + E(^j[^\mm®6)Y Using this, for any (four-marginals) 



coupling, we can rewrite (8.1 1 as 



(8.2) 



dt 



s( jf )(r))) = pE[j^'\m®m6) - Jr{»®o6) + 4''\»o®6) - Jr(««®6)], 



7(2) 



(2)^ 



so that we are dealing with four initial configurations that are aligned as follows: 



(8.3) 



•®«o 

•®o6 
•o®6 
••®6 



We perform the basic coupling, where recall the jumps are attempted at the same time at the 
same sites for the four marginals. The first and third lines contribute with a plus sign to (8.2 1 



while the others do so with a minus sign. We start with two discrepancies at the sites and 1. As 
before, under the basic coupling each of these discrepancies moves as a second-class particle, up 
until the moment when a jump is attempted from the site of one onto the site of the other. Let r 
denote this time, which could be infinity. The states of the four marginals at various times look 
as follows: 



(8.4) 



Time 

•®«6 

•®o6 
•o®6 
••®6 



Time i G (0,t) 
•a®6«c6 
•a®60c6 
•a05®c6 
•a«6®c6 



Time t— Time t 

•o®«c6 •a®«c6 •a«®c6 

•a®Oc6 •aO®c6 or •a®Oc6 

•aO®c6 •aO®c6 •a®Oc6 

•a«®c6 •a®«c6 •a«®c6 



Here a, b and c are some finite sequences of holes and first-class particles. (They are the same 
within any given column, though their values may of course differ between columns). The two 
possibilities corresponding to time r reflect the fact that the jump at time r may be from the 
discrepancy on the left onto the discrepancy on the right, or vice versa. If the jump is from left 
to right, one gets the left column; if the jump is from right to left, one gets the right column. 
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Note that before time t, the first row has one more • to the right of the ® than the second row 
does. Meanwhile the third and the fourth rows have the same number of particles to the right of 
the ®. Hence for any t < t, 



At time r, the first and the fourth rows become equal, and (since they are coupled) will still be 
equal at all later times. The same is true of the second and third rows. Hence for any t > t, 

jf ^(•©•6) - Jp'(«®o6) + Jf ^(•o®6) " jf H**®®) = 0. 

Thus we have, for any t, 

Jp^(«®«6) - Jp^(«®o6) + JpH«o®6) - JpH««®6) = i{t < t). 



t(2) 



7(2)/ 



Substituting into (8.2 1, we get 
(8.5) 



dt 



E{j[^\f,)) ^ pP{t < t) 



1 EiA^f^) 



Since {s < r} are decreasing events in s, the following limits exist: 
(8.6) 

lim P{t <t) ^ lim - P{s < r)ds = lim - / -— i;(jf )(f)))ds 

t^oo t^oo t Jq t^oo t Jq p as p t^oo t 

By Theorem |2.1[ X{ff)/t converges in distribution as i ^ oo to a Uniform random variable with 
support on [— (p — q),{p— q)] denoted again by Up. Together with the convergence to equilibrium 
in (5.2 1, this allows us to conclude that 



(8.7) 



t 



1 fip^q)~U 



in distribution. In addition, Jj (ry) is non-negative and dominated by a Poisson process of rate 
1. Hence 



lim E 

t — >oo 



-E 



^ {p-q) -u y 



Taking the limit as t ^ oo in (8.5 1 gives P(t < oo) = 1 — = as required 



3p 



□ 



Proof of Corollary |2.4| Using (8.2) and the forward equations, we have 



Pnt<r)^-E{4'\v)) 



(8.9) 



pEr,t{X{t) - 1)) - qErit{Xit) + 1), 



(2) 

since increases by one whenever the second-class particle swaps with a particle on its left, and 
decreases by one whenever the second-class particle swaps with a particle on its right. 

On the other hand, by the symmetry of the system under exchange of hole and particle and of 
left and right, rit(^X{t) — l) has the same law as 1 — rjt{X{t) -{- 1). Hence, 

(8.10) Ei^t{X{t)~\)+E7^t{X{t) + \))=\. 

As i ^ oo, the left-hand side of ( |8.9[ ) converges to pP{t — oo) which is {p — q)/3 by Theorem 
2^1 Putting and ( |8.10[ ) together then gives ([2^ and ([2j| as desired. □ 



Proof of Theorem |2.3| (b) . Consider the process started from the state •®®6. From Theorem 



2.3 'a), we know that the probability of collision of the two second-class particles is (1 +p)/3p. 



We condition on the first jump. There are three possibilities. At rate 1, the two second-class 
particles collide (either because of a jump right from site or because of a jump left from site 1). 
At rate p there is a jump right by the first-class particle at site —1, displacing the second-class 
particle at site and leading to the state •®#®6. Also at rate p there is a jump right by the 
second-class particle at site 1, leading to the state •®0®6. 
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So with probability 1/(1 + 2p) the first jump leads to immediate collision, and with probability 
2p/(l + 2p) the first jump is to one of the states •®#®6 or •®0®6. 

Notice that the dynamics of the system are invariant under the operation which reverses left 
and right and also exchanges the roles of first-class particles and holes. From this one can see 
that the probability of collision of the two second-class particles is the same from the two states 
#®#®6 and •®0®6. Let this probability be a. Conditioning on the first jump, we obtain 

l+p _ 1 2p 
3p " l + 2p ^ 1 + 2p"' 
Solving this we obtain a = (1 + 2p^)/{6p^) as required. 

9. Growth model proof 
In this section we explain the correspondence between the growth model with three types and 



the TASEP with two second-class particles, which leads to Theorem 4.1 

This is a natural extension of the correspondence between a two-type growth model and the 
TASEP with one second-class particle, developed by Ferrari and Pimentel in We start by giving 
an account of that correspondence; for a more detailed description, the reader may like to refer to 
[B] itself. 

First we make explicit the correspondence between the (one-type) corner growth model and the 
one-type tasep, as described by Rost in [TO] . 

In the corner growth model, all sites of start off empty at time 0; each site z waits until its 
two neighbours below and to the left, z — (0, 1) and z — (1, 0), have been occupied; it then waits for 
a random further amount of time which is exponential with mean 1 , and then becomes occupied 
itself. All sites {0,x) and {x,0) are already occupied at time 0. 

Write G{z) for the time that site z becomes occupied in the corner growth model. We can write 
a system of recurrences G{z) = w{z) + max{G'(z — (0, 1)), G(z — (1,0))}, for z e N"^, where we 
set the boundary conditions G{i,j) — whenever i — or j — 0. Here the 'w{z), z G are i.i.d. 
exponential random variables with mean 1. 

To represent the tasep in this way, suppose that the particles are labelled PI, P2, P3, ... 
from right to left, and the holes are labelled HI, H2, H3 from left to right, as shown in Figure |4] 
(These labels are not to be confused with priorities or with the types of particle in Section |4] the 




Figure 4. Initial condition for the tasep, with particles and holes labelled 
sequentially. The first event, at time G(l, 1), is the exchange of PI and HI. 



sequence of particles and the sequence of holes will both stay in their initial orderings throughout, 
although particles and holes will exchange places). For z = G N^, let G{i,j) be the time at 
which particle Pi overtakes hole Hj. Then exactly the same recurrences hold as in the previous 
paragraph. Once particle Pi has overtaken hole H(j — 1), and hole Kj has been overtaken by 
particle P(i — 1), particle Pi finds itself immediately on the left of hole Hj; now an exponential 
amount of time with rate 1 passes before Pi overtakes Hj. 

This gives a one-to-one correspondence between a set of states of the tasep (in which all sites 
far enough left contain particles and all sites far enough right contain holes) and a set of states of 
the growth process (with finitely many sites of occupied and satisfying the condition that if z 
is occupied then so are z — (0, 1) and z — (1,0)). G{i,j) < t means that particle Pi is to the right 
of hole Hj at time t, or correspondingly that site (i, j) is already occupied at time t. 

Ferrari and Pimentel [(? extended this correspondence to represent a tasep with a single second- 
class particle in terms of two-type growth model. 

The first element in this representation was to replace the site containing the second-class 
particle, ®, by a pair of sites, containing a hole and a particle in that order, [0#]. This pair is 
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called the *pair, consisting of the *hole and the *particle. When it is immediately to the left of a 
hole, the *particle jumps onto the hole to its right at rate 1; when this happens, we consider that 
the *pair itself moves one unit to the right, and [0#]0 becomes 0[0#]. Similarly when a particle 
on the left of the *pair jumps to the right, the *pair itself moves left, so that •[O*] becomes 
[O*]*. Hence the *pair, considered as a single unit, behaves precisely as a second-class particle. 

Hence we may represent the initial condition •®6, of a system with a single second-class 
particle, by the initial condition •[0#]6 of a one-type system. Again we label holes and particles 
sequentially, and we let G{i,j) be the first time that particle i is to the right of hole j. The 
boundary conditions are as before, except that now G{1, 1) = (since we start with particle PI 
already to the right of hole HI). 



before jump: - I 

P4 P3 P2 HI PI H2 H3 H4 

after jump: • • 



P4 P3 P2 HI H2 PI H3 H4 

Figure 5. Initial condition corresponding to the system with one second-class 
particle, and one of the two possible first jumps. The jump shown occurs at time 
G(l,2), if G(l,2) < G(2,l). 

We keep track of the labels of the *hole and *particle at time t. Let these labels be I{t) and J{t). 
We have /(O) — J(0) — 1, since initially the *pair covers hole HI and particle PI. Thereafter, 
whenever the *pair jumps right (as in Figure [s] for example), /(<) increases by one while J{t) 
remains the same; when the *pair jumps left, I{t) remains the same and J(<) increases by one. In 
particular, the position X{t) of the second-class particle at time t corresponds to I{t) — J{t) (the 
number of its jumps to the right minus the number of its jumps to the left). See Figure [5] for an 
example. 

The evolution of the process (/(i), J{t)) can be described as follows. At time t, the hole with 
label I{t) is immediately to the left of the particle with label J{t). Hence G{J{t), I{t)) < t, but 
G{J{t) + l,I{t)) > t and G{J{t),I{t) + 1) > t. The next jump in the process {I{t),J{t)) occurs 
either when particle with label J(t) jumps right, or when the hole with label I{t) moves left, 
whichever happens sooner. Let 0„ be the nth term in the jump chain of the process {I(t),J(t)). 
Then 0o = (Ij 1) and we have the following recurrence for the sequence 



(9.1) 



(0,1) if G(<^„ + (0,1)) < G(0„ + (1,O)) 
(1,0) if G((^„ + (1,0)) < G(0„ + (O,1)). 



So in terms of the growth model, behaves as follows. We start from the initial condition in 
which only site (1, 1) is occupied, and (po = (1,1). Then at each step, moves from its current 
position to whichever of its neighbours above and to the right is the first to be occupied. 

We now describe how the sequence corresponds to a competition interface in the same growth 
model when occupied sites take one of two colours. The rules are: 

• Any site {x,l), x > 1 takes the colour blue. 

• Any site {l,x), x < 1 takes the colour red. 

• For y > 1, the site (x, y) takes the same colour as whichever of its neighbours (x — 1, y) 
and (cc, y — 1) became occupied most recently. 

The colour of (1,1) is not important. 

Let B be the set of sites that take the colour blue (when they are eventually occupied) and 
let TZ be the set of set of sites that take the colour red. So TZU B — \ {(1, 1)}. Using a 



simple induction argument based on the growth rules above and on (9.1), we get the following 
characterisation which justifies the name "competition interface": the set {0„,7i > 0} consists 
precisely of those points z such that z + (0,1) £ TZ while z + (1, 0) € B. 
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In addition, if (x, y) is in the competition interface (i.e. if {x, y) = 0„ for some n, or equivalently 
if {x,y) = (/(t), J{t)) for some t) then aU the points {{x + fc,?/), fc > 0} to its right are in B and 
all the points {{x,y + k),k > 0} above it are in TZ. 

See Figure |6] for an example of the competition interface. 




X 



Figure 6. An example of the evolution of the growth model with one compe- 
tition interface. The asterisks indicate those already occupied sites which form 
part of the competition interface. For example, since the competition interface 
passes from (1,1) to (1,2), (1,3) and (2,3), we see that G'(l,2) < G(2,l), that 
G(l, 3) < G(2, 2) and that G(2, 3) < G(l, 4). 

Now we are finally ready to turn to the tasep with two second-class particles, and the cor- 
responding coexistence question in Theorem |4.1[ In the same way that the system with one 
second-class particle corresponded to a growth model with one interface, the system with two 
second-class particles will correspond to the growth model with two interfaces. 

Consider the TASEP starting from the state •®®0. We can again replace the sites containing 
second-class particles by pairs of sites, each containing a hole and a particle in that order. This 
representation works just as it did in the system with a single second-class particle, until the 
moment when one second-class particle tries to jump onto the other. At this point we move from 
[0#][0#] to 00## and we can no longer maintain the two [0#] pairs (although in fact one can 
identify a *pair in the configuration 0[0#]# and show that this *pair is in correspondence with 
a competition interface - see the remark below). In any case, up until the collision moment, the 
two second-class particle representation is valid, and in particular Theorem 2.3 'a) tells us that the 
probability that such a collision is ever attempted is 2/3. 

We will again label the particles and holes sequentially. The initial state is: 



P5 P4 P3 HI P2 H2 PI H3 H4 H5 

In this state particle PI has overtaken holes HI and H2, and particle P2 has overtaken hole HI. 
Hence we put G(l, 1) = G(2, 1) = G(l, 2) = 0. 

We will now keep track of both the *pairs: let / (t) and J (t) be respectively the labels of the 
hole and the particle covered by the left-hand *pair, and let / (t) and J (t) be those covered by 

the right-hand *pair. We have ( 7~(0), t/ (0)) (1,1) and (T(0),Tf(0)) = (2,1). More generally, 
if at time t the two * pairs are next to each other at time t if and only if 

(9.2) ~f{t) = '7{t) + 1 and t/ (i) = ~J (t) 



1. 

The local evolution both of l^I (t), J {t)j and of l^I (t), J {t)j depends on the values G{z) in just 
the same way as in the paragraph around (9.1). If we define ^ „ , n > j and ^ „ , n > j to be 
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the jump chains of these two processes, then each of them separately obeys the same recurrence 
as the process 



,) at (9.11 



Now recall the rules of the growth model with two interfaces from Section [4] In the initial 
condition, the sites (1, 1), (1,2) and (2, 1), and all other sites are unoccupied. When a new site 
becomes occupied, its colour is chosen as follows: 

• Any site (l,x) becomes dark blue. 

• Any site {x, 1) becomes light blue. 

• The site (2, 2) becomes red. 

• A site {x,y), x > 1, y > 1, {x,y) ^ (2,2), takes the colour of whichever of its parents 
became occupied more recently. 

Remark. In fact , we have elaborated slightly by introducing both dark and light blue colours at this 
point for convenience, to distinguish the clusters on either side of the red cluster. The difference 
is convenient for the argument below, but it is not significant for the result, since we only ask 
whether or not the red cluster is unbounded, and we are free to identify the two blue clusters as 
before whenever we like. In addition, once we make the distinction between the two blue clusters, 
we can consider the interface between those clusters after the red cluster disappears. This interface 
is identified by the *pair resulting from the coalescence of the *pairs related to the light-blue/red 
and red/dark-blue interfaces, if we make the jump from [0#][0#] to 0[0#]# suggested above. 

As before, the processes (^(j)nj and (^(f> nj are competition interfaces, with the following prop- 
erties. Let fe > 1 and n > 0. Then any point (/>„ + (&, 0) to the right of the lower/right interface 
is (light) blue, and any point </> „ + (0, k) which is above that interface is either red or (dark) blue. 
Similarly any point (f>n + (0; fc) which is above the upper/left interface is (dark) blue, and any 
point (j)n + (fc,0) which is to the right of that interface is either red or (light) blue. 
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Figure 7. An example of the evolution of the growth model with two interfaces. 
Here we see that G(3, 1) < G(2, 2) < G{1, 3), since the lower-right interface (indi- 
cated by asterisks) moved to (3, 1) rather than (2,2) but the upper-left interface 
(indicated by stars) moved to (2,2) rather than (1,3). In the state shown, the 
two interfaces are at (8,5) and (7,6). If G(8,6) is smaller than both G(9,5) and 
G(7, 7) then both interfaces will move next to the point (8,6). In this case the 
point (8, 6) will be red but no further points will ever become red; the red cluster 
will be surrounded. If, instead, one of (9, 5) or (7, 7) is occupied before (8, 6), the 
red cluster will grow further, perhaps unboundedly. 



From (9.2), we see that the moment (if it ever happens) when the two interfaces coincide 
for the first time is the first moment when the left-hand second-class particle attempts to jump 
onto the right-hand second-class particle. Once the interfaces meet in this way, the red cluster is 
"surrounded" and cannot grow any further. Specifically, the last paragraph tells us that all the 
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points above the upper- left interface will be (dark) blue and all those to the right of the lower-right 
interface will be (light) blue, so that no further points will be able to take on the colour red. See 
Figure [7] for an example. 

On the other hand, suppose that the second-class particles have not yet collided, so that the 
interfaces have not yet coincided. Even if the interfaces are already adjacent, it is still possible for 
them to move apart again, and any sites between the two interfaces must take the colour red; so 
certainly the red cluster is not yet surrounded. See Figure [7] for an illustrative example. 

Hence the event that the red cluster grows unboundedly corresponds to the event that the two 



second-class particles never meet. From Theorem 2.3 a), we know that that event has probability 
1/3. 
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